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ABSTRACT 

We recently reported the discovery of FORJ0332-3557, a lensed Lyman-break galaxy at z = 
3.77 in a remarkable example of strong galaxy-galaxy gravitational lensing. We present here a 
medium-resolution rest-frame UV spectrum of the source, which appears to be similar to the 
well-known Lyman-break galaxy MS1512-cB58 at z = 2.73. The spectral energy distribution 
is consistent with a stellar population of less than 30 Ma, with an extinction of Ay = 0.5 
mag and an extinction-corrected star formation rate SFRfyy of 200-300 Hjq M Q a -1 . The 
Lyman-a line exhibits a damped profile in absorption produced by a column density of about 
Nm = (2.5 ± 1.0) x 10 21 cm~ 2 , superimposed on an emission line shifted both spatially 
(0'.'5 with respect to the UV continuum source) and in velocity space (+830 km s -1 with 
respect to the low-ionisation absorption lines from its interstellar medium), a clear signature 
of outflows with an expansion velocity of about 270 kms -1 . A strong emission line from 
Hell A164.04nm indicates the presence of Wolf-Rayet stars and reinforces the interpretation 
of a very young starburst. The metallic lines indicate sub-solar abundances of elements Si, Al, 
and C in the ionised gas phase. 

Key words: galaxies: high-redshift - galaxies: starburst - galaxies: ISM - galaxies: abun- 
dances - galaxies: evolution - gravitational lensing - galaxies: stellar content 



1 INTRODUCTION on-going star formation, with strong outflows, and a dusty compo- 

nent within abundant neutral gas. 



Very few high-redshift galaxies (z > 3) have been observed spec- 
troscopically at medium and high resolution until now because they 
are generally too faint for 8/10-m class telescopes equipped with 
medium resolution spectrographs. Over the past years the most 



As of today, only a handful of galaxies are bright enough 
to be studied spectroscopically at medium and high resolution, 
and am ong them the most well-known example is MS 1512-cB58 
(cB58) jPettini et alJl200Ct 120021 ; ISavaglio et alj|2002l . z = 2.73; 

2 = 



successful technique to probe these high-redshift galaxies con- , , . x * , ■ .« , T 

. „ , ... look-back time=l 1.4 Ga). A more recent example is the 6L» .. - 

sisted in building composite spectra of low signal-to-noise sam- _ _ t , t rr- 7; ; — , t ,IUnr>cK j- j ■ n u 

1 — — B — 1 1 r 1 1 — ^-1 1 — -1 1 ; . 1 .„ 5.5 starburst Dow-Hvgelund et al. 2005) discovered in the field 

pies jShaplev et alj I2001L 12 003: Noll et al. 20041) and redshift- , , ' — ; — — T , 

r , , ' f l„ , ,IU„„J L , ,IU„„A °f a z = 1-24 galaxy cluster. The lensed galaxy FORJ0332-3557 

selected sub-samples ( Steidel et al. » 200l |;| Ando et alj | 2004|) . Very dCabanac et alj EoOS. z = 3.773, look-back time=12.1 Ga), which 

recently, high-redshift gamma-ray bursts have also led to a wealth . ^ , . ^ , ■ , , ,■ • ,■ 

is the object of the present study, is the most recent addition to this 

of low-resolution spectral data on the interstellar medium of youna , ,111,, , , , . 

1 -i f — Y select fewQ Because the galaxy is magnified 13 times by a fore- 
star-forming regions in high-redshift galaxies (Berger et al. 2009; ,11 , ,• 

— rf u »» A rr, n u»» . 1 u„„,t ^ ; — ground lens at z ~ 1, it is bright enough to be observed at medium 

Fvnbo et al. 20061: 1 Vreeswiik et alj |2004. 2006). From these com- f ,. , t , ' . °, B t , „, , ,. 

to high spectral resolution. We present here a 9h-deep medium- 



posite samples a picture emerges of very young and massive stellar 
populations of 1/1 0th- to l/4th- solar metallicities, dominated by 



resolution spectrum of the lensed galaxy taken with VFT/FORS2. 
The observations and data reduction are presented in Section 2 
while Section 3 outlines the properties of the emission and absorp- 



* Based on observations made at the ESO VLT under programs 74.A-0536 
and 78.A-0240 

f E-mail: remi.cabanac@ast.obs-mip.fr (RAC); david. vails- 1 The advent of systematic surveys on wide fie l ds is yielding many m ore 

gabaud@obspm.fr (DVG); clidman@eso.org (CL) potential candidates, see e.g. iBolton e t al. l2006); ICabanac et all txxFk) . 
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tion lines and in Section 4 we analyse the spectral energy distribu- 
tion and the underlying stellar populations. 

Throughout this paper we assume a flat FRW metric with 
Qa = 0.7, f2jvf = 0.3 and a Hubble constant normalised at 
/i70 = 70kms -1 Mpc -1 . Following the recommendations adopted 
by the IAU, the letter a is used for the non-SI unit of year. 



2 OBSERVATIONS AND DATA REDUCTION 

The observations were carried out with FORS2 on the European 
Southern Observatory VLT, under ESO programs 74.A-0536 and 
78.A-0240. A total exposure time of 8.71 hours, split into 6x 1495 
s (l"slit) and 16x1400 s (078 slit), was obtained (see Table 1). 
We used the holographic grism 600iiJ (0. 163 nm pixel -1 ) together 
with the GG435 order sorting filter, which results in resolving 
powers of 1200 and 1000 for the O'.'S the 1" slits, respectively. 
Since the O'.'S the 1" slits are slightly shifted with respect to each 
other in the focal plane, the resulting spectral coverage, 529.7- 
862.8 nm for the l"slit and 501.1-832.5 nm for the O'.'S slit, differ 
slightly. The detector was the upgraded MIT 2x4096x2048 mo- 
saic, which with the standard resolution collimator yields a physical 
scale of 0'.' 126 pixel -1 . We oriented the slits North-South and cen- 
tered them on the brightest part of the main arc (see FigureQJ. We 
used the standard 2x2 binned mode, which results in a final phys- 
ical scale of 0'.'252 pixel -1 and a sampling of 4 and 3.2 pixels for 
the 1" and 0'.'8 slits, respectively. We placed the source on CCD1 
(the upper northern chip), and used a 20" nod-on- slit strategy that 
allowed us to remove most sky artefacts. 

The science data were taken during seven nights (2 nights for 
1" slit, 5 nights for O'.'S). The ESO standard calibration scheme de- 
livers a set of bias, Hatfield and wavelength calibration (from an 
He-Ar arc) frames for each night. We, then, subtracted master bi- 
ases and divided by the 2-D normalised master flatfields. We eval- 
uated the fringing contamination to be ca. 1 % peak-to- valley in the 
most sensitive part of the frames. 

We performed two independent reductions on the 1" slit 
dataset. Firstly, we reduced separately each science exposure and 
recombined the resulting one-dimensional spectra after all calibra- 
tions. Secondly, we combined all science exposures using the FITS 
header World Coordinate info (WCS keywords) and removed the 
instrumental and sky artefacts on the combined 2-D image before 
extracting and calibrating the source spectrum. 

In the first method we used standard NOAO/IRAF routines in 
order to extract the trace with a 4th-degree Legendre polynomial 
in the dispersion direction, weighted by the spectrum variance per- 
pendicular to the dispersion direction. We removed the sky emis- 
sion from 1-D modelling of the sky emission as measured on both 
sides of the trace. We then calibrated each extracted trace with the 
wavelength calibration frames (1st order cubic spline) respecting 
the source spectrum slit position. FORS2 flexures are very small 
and the final rms of ±0.02 nm derived for the wavelength frame 
was assumed to apply to the science exposures as well. All wave- 
lengths were shifted to the heliocentric reference frame. Finally the 
combination of the resulting 6 1-D spectra (1" slit data; Table 1) 
was done using a median algorithm rejecting the minimum value 
and two maxima. 

In the second method the sky substraction was done using an 
optimal fitting while combining the 2-D spectra. The wavelength 
calibration frames proved to be stable between the two nights (1" 
slit data; Table 1). Figure|2]shows the resulting 2-D frame after sky 
substraction. 
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Figure 1. VLT/FORS2 acquisition frame (in the R band) and 1 " -slit mask 
position on FORJ0332-3557: the central object is the lens at z ~ 1 and the 
top arc is the lensed Lyman break galaxy at z ~ 3.77. 



Both methods lead to a similar final S/N and wavelength cali- 
bration accuracies across the useful wavelength range. The second 
method seems to lead to a slightly larger convolution of the ob- 
served line widths than the first method. We reduced the (/.'8 slit 
data using the second method only. 

An important calibration step is the instrument response cor- 
rection. The response curve of FORS2 was computed from spec- 
trophotometric standard observations (LTT1788 and LTT2415 us- 
ing a MOS 5" slit). Because the O'.'S and l"slits are shifted with 
regards to the MOS slit, the wavelength range of the standard cov- 
ers a slightly different range. We extrapolated the response curve 
linearly in the missing regions. 

A final step was required to remove the contribution of the el- 
liptical galaxy that produces the gravitational lens effect on the 2-D 
spectrum of the source. We used two methods, which gave similar 
results. First, we subtracted the contribution of the lens measured 
on the southern profile at a symmetrical angular distance. Second, 
we measured the dominant IR contribution of the lens over the arc 
and subtracted the best fit SED normalized to that contribution (cf. 
section[4]l. 



The ISM and stellar populations of a LBG at z=3. 8 3 




w o tn — & « i- a 



Figure 2. Continuum-normalised rest-frame ultraviolet spectrum of FORJ0332-3557 (black solid line), superimposed on its error spectrum (blue solid line, 
offset for clarity) and the spectrum of the z = 2.73 lensed galaxy cB58 (red solid line). The 2-D spectrum of the source, part of a compact lensing configuration, 
is also shown on the top. 
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Table 1. Description of spectroscopic data with VLT/F0RS2. 



Date UT time 


exp. time 


airmass 


seeing 


slit 


y-m-d h:m:s 


s 


- 


it 


n 


9004.-1 1 -1 6 01 -09-1 


1495 


1 446 


0.56 


1.0 


9004-1 1 -1 6 01 -98 01 


1495 


1 326 


0.55 


1.0 


9004-1 1 -1 6 01-54-99 


1495 


1 231 


65 


1 


9004 11 16 09-90 00 


1495 


1 161 


lf\ 
u. tv 


1 

1 At 


900S-01 -^0 01 -90-4 S 


1495 


1 103 


50 


1.0 


900S 01 10 01-4.6-96 


1495 


1 154 


54 


1 

l At 


9006-1 0-16 06 09-90 


1400 


1 028 


0.66 


8 


9006 10 16 06-96-S7 


1400 


1 090 


53 


o ^ 

U.o 


9006-1 0-1 6 06-^1 -XK 


1400 


1 021 


0.61 


0.8 


2006-10-16 07:15:52 


1400 


1.031 


0.50 


0.8 


2006-10-16 07:55:52 


1400 


1.068 


0.73 


0.8 


2006-10-16 08:20:31 


1400 


1.104 


0.61 


0.8 


2006-10-21 07:37:59 


1400 


1.070 


0.48 


0.8 


2006-10-21 08:02:09 


1400 


1.106 


0.61 


0.8 


2006-11-21 07:03:26 


1400 


1.266 


0.85 


0.8 


2006-11-22 04:33:11 


1400 


1.023 


0.85 


0.8 


2006-11-22 04:57:20 


1400 


1.036 


0.89 


0.8 


2006-11-28 04:24:07 


1400 


1.030 


0.54 


0.8 


2006-11-28 04:48:36 


1400 


1.048 


0.63 


0.8 


2006-11-28 05:20:59 


1400 


1.087 


0.50 


0.8 


2006-11-28 05:45:14 


1400 


1.130 


0.57 


0.8 


2007-01-25 04:30:28 


1400 


1.916 


0.67 


0.8 



3 EMISSION AND ABSORPTION LINES 

3.1 Source redshift 

In the low-resolution discovery spectrum (Cabanac et al. . 2 005T) we 
estimated a source redshift of z = 3.773 ± 0.003 from a few 
strong interstellar lines. The improved signal-to-noise obtained 
here (continuum peak S/N~15) allows us to measure here a re- 
vised systemic redshift of zlis = 3.7723 ± 0.0005, using 11 
stronger detections of 5 atomic species : C, O, Al, Si, and Fe, 
and excluding uncertain identifications and blended lines (see Ta- 
ble 2). We marginally detect two important photospheric lines 
O IV A 134.33) and S V A 150.18 (see Fig. Hi. Hence we can only 
make a very tentative estimate of a possible velocity offset between 
the photospheric lines and the low-ionisation ISM lines of about 
+ 110 ± 30 k m s~ . This is about half t he value inferred in cB58 of 
+210 km s^ 1 JPettini et alj2000l.l2002h . 

3.2 The Lyman-a lines 

Figure [2] shows the 2-D spectrum and the associated 1-D spectrum 
normalised to the continuum. Bluewards of Lya the spectrum is 
poorly defined mostly because of foreground hydrogen absorption, 
and the continuum has been assumed to be located at the peaks of 
the spectrum. We also added the low-res olution VLT/F ORS2 spec- 
trum at wavelengths shorter than 110 nm dCabanac et al.l2005h . The 
Lya, Ly/3, and Ly7 absorption features are conspicuous, and unre- 
solved Ly<5 absorption is detected as well. 

Figure [3] zooms on the region around Lya. A faint emission 
peak at Av — +830km/s of the Lya absorption centroid and 
~ 2 ± 0.5 pixels (C//51 ± CK'13) north of the continuum emis- 
sion is clearly visible. This feature can be interpreted either as a 
very bright Hll region in the periphery of the main galaxy, or, more 
probably, as the signature of an expanding outflow. Assuming that 
the results of detailed 3 -D Monte Carlo Lya radiative transfer codes 
IVerhammee t al. ( 2006) (where the Lya emission is produced by 




Figure 3. Zoom on the 2-D spectrum around Lya, showing the spatial offset 
of ~ 2 ± 0.5 pixels (0.51 ± 0.13arcsec) and the velocity shift of Av = 
+830 kms^ 1 , indicative of outflows (121.567 nm position is the narrow 
green vertical line). The lower trace is the spectrum of the elliptical galaxy 
that produces the gravitational leasing effect. 



an expanding shell that is being ionised) can be applied to this case, 
we can infer that the expansion velocity of the shell that gives rise 
to the Lya emission is constrained to be between V exp ~ Av ~ 
800 km s -1 (low column density case) and V exp ~ An/3 ~ 270 
km s~ 1 (large column density case). It is therefore essential to mea- 
sure the Hi column density in this context. 

If we assume that the Lya absorption can be fitted by a 
damped Voigt profilfl we derive a column density of Nm = 
(2.5 + 1.0) x 10 21 cm -2 . Figure|4]shows the best-fit damped Voigt 
profile and errors, along with the emission profile in velocity space 
obtained when the absorption profile is subtracted. The column 
densit y is remarkably large, about 4 times larger than in MS1215- 
cB58 jPettini et alj |2002| - ISavaglio et all |2002|) . If we adopt the 
mode l of an expanding shell at constant velocity dVerhamme et all 
2006), this column density still implies a large expansion veloc- 
ity of about V exv ^ 270 km /s, similar to the one found in cB58 
( Schaerer & Verhammell2008 |) , and at the upper end of the Lyman- 
break galaxies dShaplev et ai]|2003l) if they do have large Hi col- 
umn densities. It is remarkable that in spite of the very large column 
densities, such strong Lyman-a emissions are detected, confirming 
that the interstellar medium is very clumpy, yielding smaller ef- 
fective optical depths which a re very similar to th e ones inferred 
in the low-redshift starbursts dVairs-Gabau"3ll993l) . and allowing 
many Lya photons to escape. 

The neutral gas content derived with Voigt profile fitting is 
fully consistent with the neutral gas content derived independently 
from the analysis of stellar populations (Section|4j with the gas-to- 
dust ratio taking Milky Way or LMC values, but is ~ 5-10 times 
smaller than the values observed in the SMC. 



2 This is clearly a rough approximation since the spectrum comes from the 
integrated light of the source and there are photons that are scattered into 
the line of sight. Hence the profile cannot be a pure absorption profile, and 
the inferred column density is just indicative. 
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Figure 4. The best-fit Voigt profile corresponding to a column density of 
(2.5 ± 1) X 10 21 cm -2 (solid and red dotted lines) is shown on the rest- 
frame spectrum (top frame). The bottom frames show the same fit in ve- 
locity space Av. On the panel at the bottom right the model is subtracted 
out of the spectrum, showing a conspicuous emission feature peaking at 
~ 830 km/s. This feature is also offset spatially by 2 pixels (ca. 0'.'5) from 
the main UV-continuum emitting source (see Fig.[5J. 



3.3 The Helium II A164.04nm emission line: the signature of 
very massive stars 

A broad emission feature at 164.04 nm is seen in Fig. [2] at a red- 
shift of z = 3.7738 consistent with other faint photospheric fea- 
tures. We associate this feature with the well-known Hell emis- 
sion produced by Wolf-Rayet stars. Although not present in cB58, 
this feature has alrea dy been detected in absorption in young LBGs 
dShaplev et af]|200lL ; cf. Fi g. 15 and 16) and in emission in the 
composite LGB spectrum of Sh aplev et al.l J2003I) but the origin of 
the feature cannot be asserted by the authors because the compos- 
ite spectrum mixes a wide range of stellar populations and ages. If 
FORJ0332-3557 is similar to other LBGs, it belongs to the strongly 
absorbed Lya quartile, and the case for a very young stellar popula- 
tion (~5 Ma) of massive stars is strong. The He II feature of Fig. [2] 
shows the typical asymmetric profile produced by strong stellar 
outflows of the evolved descendants of O stars more massive than 
M > 2O-3OM . The equi valent width Wo (Hell) y ields the ratio 
of WR to O stars following Sch aerer & Vaccal dl998 



log 



WR 



WR+O 



which for Wo (Hell) = 0.22 ± 0.02nm implies a ratio 
OfWR— 5 ± 1. In turn, this value sets a strong age limit on 
the ste llar population of ^6 Ma ra t her independent ly of metal - 
licity dCervino & Mas-Hessd 1 19941 : ISchaeren 120031) . Stochastic 
effects due to the sampling o f the stellar mass fu nction (e.g. 
ICervino & Valls-Gabaudl d2003h : ICerviSo et ail J2003h ) imply that 
only the most recent burst can be dated through this short-lived fea- 
ture. 



3.4 Interstellar absorption lines 

The spectrum of FORJ0332-3557 shows characteristic absorption 
lines from starburst galaxies. They are summar ised in Table 2. The 
medium-resolution spectrum of the LBG cB58 |Pettini et al.[ ( 2000) 
is shown in Fig.|2]in red, superimposed on our source, for a direct 



comparison of the insterstellar features. It is immediately obvious 
that the FORJ0332-3557 source is qualitatively similar to cB58. 

All common interstellar absorptions are found in FORJ0332- 
3557. The strong absorption features include low-ionisation lines 
associated with neutral gas (Sill AAA126.04 130.47 152.67, ClI 
A133.45, Ol A130.22, Alii A167.08, Fell A160.84) and high- 
ionisation lines associated with a hot gas phase (Si IV AA139.38 
140.28, CIVAA154.82 155.08). Table 2 lists the ion line iden- 
tification, vacuum rest-frame wavelength A^ b c , observed wave- 
length \° b a B c , redshift z, rest-frame equivalent width Wo, oscilla- 
tor strengths /, column density, ion abundance with respect to so- 
lar [X/H]q and comments. Additional uncertain identifications 
are question-marked, the lines noted i? might belong to interlop- 
ers at unknown redshift(s). We emphasize that the derived Wo are 
very sensitive to both sky subtraction and continuum normalisation, 
hence the systematic errors caused by the continuum normalisation 
have tentatively been computed and are shown to be close to pho- 
ton counting errors, while the sky subtraction error is much more 
difficult to quantify. As a sanity check, we computed the equiva- 
lent widths of absorption lines in the spectrum of cB58 and found 
our measurements to be fully consistent with those published by 
IPettini et ail fc002t) . 

No nebular emission lines are detected in the present spectrum 
other than Hell A164.04nm (§ !3.3t . A weak detection o f Cm] A 
190.87 seen on a previous spectrum JCabanacetaIB2 005) suggests 
that contamination of the high-ionisation lines by nebular emission 
is present but small. P-Cygni profiles are visible on C II A133.45, 
andCivA155.08. 

There are several ways to de rive the abundances in the inter- 
stellar medium of dis tant galaxies dSpitzeJl978MPettini et alj2002t 
ISavaglio etafll2002l) . Ideally one should build a curve of growth 
by fitting Voigt profiles and Doppler parameters b for all ions in- 
dependently. Because the resolution of our observed spectrum is 
just under the resolution one needs for Voigt profile fitting, and is 
penalised by a low signal-to-noise ratio, most of the strong lines 
appear saturated, and most weak lines are dominated by noise. 

A careful analysis of the ISM metallicity goes beyond the 
present paper and will be done elsewhere. Here we present only 
qualitative arguments on the curve of growth, and Doppler param- 
eters b. Assuming that the interstellar medium in FORJ0332-3557 
is optically thin, one can infer lower limits to column densities, 
log(iV[cm -2 ]), and abundances (given in Table 2) by taking the 
optically thin approximation 



= (-2.83±0.03)+(1.52±0.05) log[W (Heii)] ,(1) logiv[cm - 2 ] 



19.053 + log 



Wa 



A 2 / 



(2) 



where / is the line oscillator strength. The equivalent width, W\, 
and the wavelength, A, are in nm. A tentative curve of growth 
indicates that the ion abundances could be 2-3 dex larger for a 
Doppler parameter of b = 50kms _1 . In this context, the most 
constraining line, besides Sill* A153.3 which may be blended, is 
FeIlA160.8, which appears unsaturated and whose small equiva- 
lent width is similar to the one measured in cB58 and would yield 
b ~ GOkms -1 , similar to the b ~ 70 km reported in cB58 
jPettini et all2002l) . 

Compared to cB58, FORJ0332-3557 W are lower by factors 
of 2-3 (CIVAA155.08 155.08, AlnA167.08, O IA130.22) to a fac- 
tor of 1-1.2 (SHVAA139.38 140.28, FeIlA160.84). More detailed 
analyses on elemental abundances and depletion in the interstellar 
medium of FORJ0332-3557 will be presented elsewhere. 



6 R.A. Cabanac, D. Valls-Gabaud and C. Lidman 



Table 2. Interstellar absorption lines 



Ion 


\ lab 
"vac 

nm 


\obs 

vac 

nm 


redshift 


nm 


/ 


logTvi 
cm -2 






Comments 


n i L.y f 








l.OV IIZ U.OU Q JQ 


u.uzyu 










rl l L*yp 


1 09 






2.1U ± U.»U_ 6Q 


fi 07Q 1 9 

U.U iy 












1215 7 








4164 


21.4 ±0.2 






He II 


164.04 


783.09 


3.7738 


u.zz =t u.ui_ Q2 










WR fpatnrp 

VV IV ill III \. 


C II 


133.45 


637.01 


3. 7732 


21 + rn+ ' 03 


0.1278 


15.019 ±0.10 


> —2.8 


blended with C II*A133.57 


C iv 


154.82 


738.74 


3.7716 


nin + n 00+006 


0.1908 


14.40 ± 


21 


> —3.6 




C iv 


155.08 


740.00 


3. 7718 


19 + 05 +0 ' 02 


0.09522 


14.97 ±0 


25 


> —3.0 


sky contamination 


Nl? 


~ 120 


572.70 


3. 7718 


Dl + n 10+0. os 


0.04023 


15.84 ±0 


25 


— 1.6? 


triplet NlA 119.95 120.02 120.07 


Nm? 


132.43 


632.24 


3. 774 1 












blended with CIIA132.39 


Ol 


130.22 


621.31 


3. 7714 


0.15 ± 0.03+225 

— 0.03 


0.04887 


15.30 ±0 


17 


> -3.0 


blended with Si II A 130.44 


Oiv 


134.34 


641.27 


3.7735 


0.14 ± 0.01 + n '„5 
— 0.03 










photospheric 


Alii 


167.08 


797.41 


3.7726 


0.19 ± 0.02+J1!!! 

—0.03 


1.833 


13.62 ±0 


22 


> -2.5 


sky contamination 


Sl/i? 


138.16 


659.29 


3.7719 


0.11 ±0.01+°'°2 












S Il/i? 


125.38 


598.57 


3.7740 


0.12±0.02±g;°* 


0.01088 


15.90 ±0 


15 


-1.1? 




S Il/i? 


125.95 


601.31 


3.7742 


0.19 ± 0.04+g'° 3 


0.01624 








blended with Sin A126.04 


Sv? 


150.18 


716.39 


3.7741 


0.3 ± 0.03t° ° 4 


0.00545 








photospheric 


S Il/i? 


151.12 


721.43 


3.7740 


0.07±0.0i1q: 2 










blended with SinA151.21 


Sill 


126.04 


601.84 


3.7749 


0.19 ± 0.04±°;° 3 


1.007 


14.13 ±0 


11 


> -2.8 


blended with S IIA125.95 


Sill 


130.47 


622.49 


3.7710 


0.16±0.0ll°;°^ 


0.094 


15.05 ±0 


37 


-1.9? 


blended with O IA130.22 


Sill 


152.67 


728.39 


3.7710 


0.24±0.03±°; 2 


0.130 


14.96 ±0 


11 


> -2.0 




Sin*? 


153.34 


731.40 


3.7697 


0.04±0.01+°;°j 


0.132 


14.20 ±0 


10 


-2.8? 




Si IV 


139.38 


665.09 


3.7718 


0.19 ± 0.03l°;p2 


0.5140 


14.33 ±0 


10 


> -2.6 




Si IV 


140.28 


669.38 


3.7718 


0.21 ± 0.03±°;°| 


0.2553 


14.66 ±0 


12 


> -2.3 


blended with S IA140.15? 


Fell 


160.84 


767.68 


3.7728 


0.12 ± 0.02to; 02 


0.058 


14.95 ±0 


11 


-2.0 





a Errors are given as statistical (±photon noise) and systematic (_ systematic uncertainty in location of the continuum) 
^ Lower limits based on the assum ption of an optically th in medium 
c Assuming solar abundances from Asplund et al. 1 2 0051) 



4 STELLAR POPULATIONS 

After analysing the absorption and emission lines here we study the 
properties of the continuum, using recent stellar synthesis codes, 
which, for the first time, predict high-resolution spectra of star- 
burst populations. We used both STARBURST99 ( Leither er et all 
19991: IVazauez & Leifhere3 120051) and sed@ dBuzzonil 1 19891 : 
Cervino et al. 20021 : Cervino & Lurid mn3l2004h n an attempt to 
derive plausible combinations of ages, metallicities and extinctions 
by dust from the FORJ0332-3557 ultraviolet spectrum (Fig.f2j- 



4.1 Spectra predicted by starburst99 

We build two series of STARBURST99 spectra using (1) an IMF 
with a Salpeter slope from 1-100 Mg0, (2) a self-consistent asso- 
ciation of low metallicity tracks with LMC/SMC atmospheric stel- 
lar templates, and (3) solar metallicity tracks with galactic atmo- 
spheric stellar templates, from ages 0.01 Ma to 50 Ma, for metal- 
licities Z of 0.001 (1/2OZ ), 0.004 (1/5Z ), 0.008 (2/5Z ), and 
0.020 (Zq). All STARBURST99 models were fitted to the spec- 
trum of Fig. [2] pre-filtered with a Savitsky-Golay filter (conserv- 
ing the equivalent-width, FWHM and line profiles), to the same 
resolution as the STARBURST99 synthetic templates. The region 



3 Tests where the upper stellar mass cutoff of the IMF is is reduced to 
30 Mq give very similar best-fit ages. 



below 123 nm was excluded from the fits, because the source is 
heavily contaminated by foreground interlopers. For solar metallic- 
ities, STARBURST99 high-resolution atmospheric templates go up 
to 180nm, but only 160nm for subsolar SMC/LMC metallicities. 
Only the highly ionised absorption lines were considered for the 
fitting, to avoid excessive interstellar contamination. An additional 
range of wavelengths (near 137 and 144 nm) was excluded due to 
telluric contamination. The regions included in the fit are shown 
as solid lines in Fig. [5] Different sets of constraints can be derived 
by fitting the line equivalent widths of the normalized spectrum and 
by fitting the shape of the spectrum continuum. Fig.|5]compares the 
best fit models with the normalised spectrum for two extreme sce- 
narios of star formation: (i) a single burst of star formation (bottom) 
and (ii) a continuous SFR (top). The best fit models for Z=0.004 
have ages of 8 Ma (single burst) and 29 Ma (constant SFR). The 
age of the best fit models is found to be insensitive to metallicity 
in the case of single bursts. In the case of a constant SFR, increas- 
ing the metallicity tends to increase the age of the best fit model : 
29 Ma at .2=0.001 , > 49 Ma at Z=0.008, and more than 50 Ma at 
Z=0.02. We explored solar metallicity models with ages of 200 Ma 
and although the reduced x" 1 decreased slightly the highly-ionised 
absorption lines did not get deeper significantly. The reduced \ 2 
shows that constant SFR models become insensitive to age above 
20 Ma, whereas single burst models appear to isolate a well-defined 
small region (6-9 Ma). However, neither of these scenarios can re- 
produce the observed depth of the absorption lines, which probably 
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argues for interstellar origin for the highly ionised lines. Finally, 
the single burst scenario seems to produce slightly deeper absorp- 
tion lines, whereas the constant SFR scenario better fits the P-Cygni 
profiles of CIVAA155.08 155.08 nm and the red side of Lya. We 
also tested different stellar wind models without any significant dif- 
ferences. 



4.2 Spectra predicted by sed @ 

The stellar population synthesis code SED@ is included in the 
Legacy Tool project of the Violent Star Formation Netwo rk 
teuzzonilll98"9l : ICervino etai]|2002l ; lc"ervino & Luridianall2004) . It 
computes UV spectra with the highest resolution currently avail- 
able, consistently computing the errors introduced by a limited 
number of bright stars in a small burst (this effect can probably 
be neglected in our case unless the lens amplifies a single star- 
forming region within the source). We use a series of instantaneous 
burst models from to 10 Ga, a Salpeter IMF, and Geneva stel- 
lar tracks at 4 metallicities: Z of 0.001 (1/20Z Q ), 0.004 (l/5Z Q ), 
0.008 (2/5Z®), and 0.020 (Zq). As in the previous section, we treat 
independently the constraints coming from the normalised UV ab- 
sorption features, and the overall photometric SED of the source. 
For the UV-spectrum fitting, we isolate the strongly-ionised absorp- 
tion features containing most of the signal and exclude the weakly- 
ionised features dominated by interstellar features, using the wave- 
length intervals 120-125, 127-129, 131-133, 134-143, 145-152, 
and 153-158 nm. We normalise each interval to the continuum by 
dividing each wavelength bin by the linear interpolation between 
the two extreme values of the interval. We compute the reduced x 2 
between the templates and the source UV spectrum for a grid of ex- 
tinctions (0 < Ay < 1), and ages (1-1000 Ma). The photometric 
SED of the templates are integrated over the instrumental transmis- 
sion curves (filter+telescope+CCD detector) given by VLT/FORS 1 
(B, V and R c bands) and VLT/ISAAC (J, H, K s bands) man- 
uals, normalised to VEG A magnitude using Kurucz's solar model 
(Castelli & Kurucz 1994). Figure[6]shows the resulting reduced \ 2 
contours. Ages are along the horizontal axis (steps of 1 Ma from 

to 10 Ma, steps of 5 Ma from 10 to 100, steps of 100 Ma from 
100 Ma and older ages). The extinction Ay is given on the verti- 
cal axis. The left panel shows the contours of \ 2 values obtained 
in the fitting of the UV features described above (the contour lev- 
els are given at 1, 3 and 5% above the value of the minimum \ 2 )- 
The central panel shows the contours obtained in the fitting of the 
photometric SED (contour levels at 1.5, 3, and 5 times the min- 
imum x 2 ). while the right panel shows the combined constraints 
from both the overall SED and the UV features (contour levels at 

1 3, and 10% above the minimum x 2 )- The left frame shows a 
complete degeneracy in extinction, as expected since we normalise 
the continuum and the UV is no longer sensitive to extinction. The 
central panel shows the well-known degeneracy between age and 
extinction but clearly shows a tendency to prefer young ages with 
higher extinction, rather than older ages with small extinctions. It 
should be noted that the \ 2 °f the UV fitting alone shows a mini- 
mum at Z=0.004, whereas the y 2 of SED-fitting alone favours the 
smallest metallicity models Z=Q.QQ\. The combined set leads to 
Z=0.008 (Fig.[6]l but it is clearly not a very constraining result due 



4 See the SED@ online Reference Manual at 
|http : / /www . iaa ■ es/ ~mcs/sed@ /index .html| for more in- 
formation. 



the degeneracies involvecjf]. Because the two sets of constraints are 
independent and almost orthogonal, the combined \ 2 shown in the 
right panel of Fig. [6] (each \ 2 is given the same weight) gives a 
tighter constraint on the age,- Ay parameter space. The combined 
best fit parameters (Z, age [Ma], Av[mag]) are (0.001, 15, 0.62), 
(0.004, 20, 0.54), (0.008, 20, 0.54), (0.020, 10, 0.64), (0.040, 15, 
0.52) for reduced \ 2 values of 0.853, 0.800, 0.788, 0.792, 0.789. 

In conclusion, both stellar synthesis codes give a consistent 
picture of the global properties of the source galaxy in FORJ0332- 
3557, although they marginally disagree on details. They both find 
young stellar populations of < 30 Ma and the presence of some 
dust extinction. Neither the metallicity, nor the IMF can be robustly 
constrained although SED@ marginally favours Z=0.008. We em- 
phasize that none of the models, at any metallicity and age, can fit 
the depth of the observed strongly-ionized absorption lines. STAR- 
BURST99 points to a young stellar population of 8-10 Ma, in the 
case of an instantaneous burst, and a few tens of Ma, in case of a 
constant SFR, whereas SED@ models seem to favour slightly older 
ages (10-20 Ma) for an instantaneous burst. Lastly, STARBURST99 
tends to give a better fit of the absorption features than SED@, in- 
cluding interstellar features of course because SED @ uses synthetic 
stellar spectra whereas STARBURST99 uses observed stellar tem- 
plates where the Milky Way interstellar absorption features have 
not been removed. Finally, SED @ over predicts an absorption line 
Si II A 126.5 nm and a P-Cygni profile of OlA 130.22 nm. Finally, 
the ages derived by stellar synthesis are older than the presence of 
the Helium II emission line seems to suggest. This is not contradic- 
tory but rather suggests that the two analyses are not sensitive to the 
same features, which probably come from different regions of the 
galaxy and may point at a range of ages and multiple star formation 
episodes. 

4.3 Integrated spectral energy distribution and dust 
extinction 

The photometry in optical bands B B e Sse i, Vessel, flcousin, ^Guim, 
and in i nfrared bands Js , H , and K s is described in Caba nac et al.l 
(2005): IValls-Gabaud et alj (2006). Because the deflector is com- 
pact in the optical range, the source i?cousin and iGunn magnitudes 
were simply extracted in an annulus of inner/outer radii of l'.'5/3". 
The source is faint in BbcsscI and Vsessei bands, and we can only 
put an upper limit to its flux in these bands. The IR magnitudes 
are more difficult to extract because the foreground deflector cov- 
ers part of the sou rce. We subtracted the central len s with GALFIT 
dPeng et alj|2002l) . using the best fit parameters of Caba nac et al.l 
(2005). The derived IR magnitudes are sensitive to the lens sub- 
straction and we estimate that the systematic errors are of the order 
of 0.3 mag. We obtain the integrated Spectral Energy Distribution 
(SED, in Vega magnitudes) as shown in Figure|7j 

We derived dust exti nctions Ay, and ages o f the stellar pop- 
ulations using HyperZ teolzonella et alj bOQO). for the known 
redshift of 3.7723, an extinction law similar to the Small Mag- 
ellanic Cloud and several libraries of HyperZ templates, span- 
ning single burst to constant SFR. The stellar Initial Mass Function 
(IMF) is that of iMiller & Scald d 19791) and the metallicity of Hy- 
perZ templates is fixed at the solar value. Although it is likely that 
FORI0332-3557 has a subsolar metallicity at a look-back time of 
~ 12 Ga, at face value solar metallicities are as good a guess as 

5 See for example Schaerer & Pelld ( 2005) for other examples of degen- 
eracies in the SED fitting of high-z galaxies. 
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SFR=Constant Z ,=0.004 IMF=2.35 Age = 29Ma 

met & 



120 125 130 135 140 145 150 155 160 

wavelength [nm] 




SFR=Instant Z ,=0.004 IMF=2.35 Age = 8.01Ma 

met fa 



120 125 130 135 140 145 150 155 160 

wavelength [nm] 

Figure 5. The best fit STARBURST99 models (red lines) are overlaid on the normalised spectrum of FORJ0332-3557 source (blue line) for two extreme 
scenarios of star formations with a Salpeter IMF from 1-100 Mq, at a metallicity Z=0.004 (1/5 Zq). The parts of the spectrum excluded of the fit are shown 
as dotted lines (see text). The top frame shows a constant SFR for 29 Ma. Ages older than 20 Ma are strongly favoured. The bottom frame shows an 8 Ma-old 
instantaneous burst of star formation. No model (at any available metallicities (1/20-1 Zq) can reproduce the observed depths of the lines. For a constant SFR, 
larger metallicities produce older ages (29 Ma at Z = 0.001, >49Ma at Z = 0.008). Instantaneous burst models tend to produce deeper lines while constant 
SFR models tend to better fit the observed P-Cygni profiles of C IVA155.08 nm. 




5 10 35 60 85 300 5 10 35 60 85 300 5 10 35 60 85 300 

Age [Ma] 



Figure 6. Reduced \ 2 contours are given for the best fit Z=0.008 SED@ models for a grid of ages and extinctions. The contour levels are arbitrarily chosen 
to be 1, 3 and 10% of the reduced \ 2 minimum, to outline the trends. The left panel shows the contours derived for a continuum-normalised set of deep 
absorption lines alone (cf. text), the center panel shows the constraints derived from the photometric colours alone, and the right panel shows the combined 
contours of the two independent sets taken together. The optimal parameter set is age=20±5 Ma, Ay = 0.55 ± 0.02. 



sub-solar metallicities. In order to estimate the effect of metallic- 
ity on A v , we also fed HyperZ with the best fit STARBURST99 
templates (cf. previous section) of Z = 0.004. 

Both Single-burst and Irregular HyperZ templates yield the 
same dust extinction Av = 0.5 mag (for Rv = 3, E[B — V] = 
0.17) for the best fit 23 Ma-old single-burst template and the best 
fit 128 Ma constant SFR template (Figure[7]l. The ages of the stel- 
lar populations are not well constrained by stellar synthesis, and 
Fig. [7] illustrates this point. One template is 100 Ma older than 
the other and yet they are indistinguishable with broadband data 
only. The inferred col our excess of E(B — V) = 0.17 is similar 
to the value found bv lShaplev etafl J2003I) (0.169 ± 0.006) for the 



Group 1 (Lya deficient) LGBs, using a similar extinction law and 
photometric bands, but a 300-Ma-old constant star formation rate. 
FORJ0332-3557 also falls in t he middle of the ex tinction distribu- 
tion of the z ~ 4 LBG sample Isteidel et alj Jl999l> . 

Following Petti ni et"aT] ( feOOOh . we infer the extinction at 
150 nm, Also and the Hi content, N(Hl ). For an SMC like extinc- 
tion law, A150 = 12.6 x E(B — V) ~ 2 - 2.5. Assuming gas-to- 
dust ratios of< N{m)/E{B - V) >= 4.93 x 10 21 cm~ 2 mag _1 
for the Milky Way, 1 - 2 x 10 22 cm~ 2 mag _1 for LMC, and 
5 x 10 22 cm- 2 mag~ 1 for the SMC, we find 7.6 x 10 20 cm~ 2 < 
AT (Hi ) < 10 21 cm" 2 (MW), 1.5 x 10 21 cm~ 2 < AT (Hi ) < 4 x 
10 21 cm -2 (LMC), and 7.5 x 10 21 cm~ 2 < AT (Hi ) < 10 22 cm~ 2 
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Figure 7. The BVRIJHK S SED of FORJ0332-3557 (red crosses, the 
vertical bars indicating systematic errors and the horizontal ones giving the 
spectral bandwith of each filter). The top frame shows HYPERZ best fit solar 
metallicity templates: single-burst template of age 23 Ma (blue solid line) 
and constant SFR followed by passive evolution template of age 128 Ma 
(green dotted line). Both templates are the best fit with an extinction Ay = 
0.5 mag. The bottom frame shows HYPERZ best fit with STARBURST99 
Z=0.004 templates: instantaneous burst of age 8 Ma (blue solid line) and 
constant SFR of age 35 Ma (green dotted line), both with a best fit extinction 
of Ay = 0.6 mag. The bottom frame shows a clear difference between a 
constant SFR and an instantaneous burst near the 400 nm break. We can 
favour an instantaneous burst scenario based on Kg band data. 



(SMC). The column density inferred from the Lya absorption line 
is compatible with the first two values and would be inconsistent 
with an SMC gas-to-dust ratio. 



4.4 Dust-corrected star formation rate 

The UV (150 nm) luminosity is a standard indicator of (massive) 
star formation. Witho ut correcting for dust extinction and assum- 
ing a constant SFR, ICabanac et alj d2005h derived SFRuv ~ 
31 (A/12.9) -1 h^ M©a _1 , where A is the gravitational ampli- 
fication produced by the lens. Using the extinction measured in 
the previous section, A150 = 2 — 2.5, we get SFRuv ~ 196 — 
310 (A/12.9)" 1 M©a _1 . Using the same method on cB58, 
IPettini et all l l2000h derives SFRuv = 395 M a -1 (for H = 
70 km s" 1 Mpc" 1 , 50 = 0.1). This puts FORJ0332-3557 in a very 
dynamic star forming episode, s imilar to what is observed in SDSS 
Jl 147-0250 teentz et aljiool) . 



5 CONCLUSIONS 

We have presented a medium-resolution spectrum of FORJ0332- 
3557, a lensed z = 3.7723 ± 0.0005 starburst galaxy, magnified 
13 times, very similar to both z = 2.73 cB58 and the composite 
spectrum of Lyman-break galaxies at z ~ 3. 

We find spectral signatures of outflows of ~ 270 km s^ 1 , 
commonly found in starburst galaxies, through a distinct residual 
Lya emission off-centered by 0'.'5 and peaking at ca. +830 km s" 1 . 
This offset emission is most probably produced by the outflow 
rather than by an isolated, bright Hll region within the galaxy. 

By combining visible-to-infrared colours and spectral fea- 
tures, we derive consistent ages of S3 10 Ma for single bursts and 
20-40 Ma for constant SFRs and an extinction of Ay = 0.5 us- 
ing two independent stellar populations synthesis codes. Young 
stellar ages are also supported by the detection of the Wolf-Rayet 
Hell A 1 64.04 nm line, indicator of stellar ages in the range 4-6 Ma. 
There is a clear trend to sub-solar metallicities in both cases. The 
ongoing rate of star formation, computed from the 150 nm contin- 
uum is SFRuv « 200 - 300 M© a" 1 . 

We also derive preliminary lower limits to the abundances 
of some the low-ionisation interstellar lines suggesting sub-solar 
metallicities in the ionised gas phase. 
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